Crossover from Conserving to Lossy Transport in Circular Random Matrix Ensembles 
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In a quantum dot with three leads the transmission matrix ti2 between two of these leads is a 
truncation of a unitary scattering matrix S, which we treat as random. As the number of channels 
in the third lead is increased, the constraints from the symmetry of S become less stringent and ti2 
becomes closer to a matrix of complex Gaussian random numbers with no constraints. We consider 
the distribution of the singular values of ti2, which is related to a number of physical quantities. 



Random matrix theory has enjoyed enormous success 
in a wide range of fields of physicsjj. One impressive 
success has been in mesoscopic physics [3, where the scat- 
tering matrix between two leads (a source and a sink) of 
a disordered quantum dot has been modelled as a ran- 
dom matrix. Such scattering matrices can be from the 
circular unitary ensemble (CUE or unitary matrices) if 
time reversal symmetry is broken, the circular orthogonal 
ensemble (COE or unitary symmetric matrices) if time 
reversal is unbroken, or the circular symplectic ensemble 
(CSE or unitary self dual matrices) if spin-orbit scat- 
tering is present and time reversal is unbroken^. An- 
other success of random matrix theory has been in the 
arena of radio communicationsi where transmission 
matrices between transmitting and receiving antennas 
have been modelled as random matrices with all matrix 
elements being independent and identically distributed 
(i.i.d.) complex Gaussian. These two examples are, in 
fact, quite closely related to each other in that both stem 
from wave transmission problems. The fundamental dif- 
ference between these two cases is that the former is "con- 
serving" while the latter is "lossy" . When an electron 
enters a dot, it must exit one of the two leads whereas in 
the radio communication problem most of the photons 
which are transmitted are lost to the environment and 
neither reach the receiver nor return to the transmitter. 
In this manuscript we will study the crossover between 
these two cases by studying a quantum dot with a "third 
lead" grounded at the same voltage as the sink of the dot. 
As the number of channels to ground through the third 
lead is increased (i.e., the amount of loss is increased), 
the statistics of transmission between the source and the 
sink become more and more similar to the i.i.d. complex 
Gaussian case independent of the underlying symmetry 
(GOE, CUE, or CSE). We note that this application of a 
"third lead" is quite distinct from that used previously in 
the mesoscopic literature to model dephasing|2]. In that 
case, a constraint is imposed that no current may flow 
out the third lead; whereas in the present case, we allow 
current to flow out the third lead, giving us true loss of 
electrons. We will not consider dephasing here, so this 
work strictly applies to dots at zero temperature. 

We thus consider a quantum dot with three leads, a 
source (lead 1) at voltage V, a grounded sink (lead 2) 



at voltage 0, and a loss lead (lead 3) also grounded at 
voltage 0. We assume there are Ni fully open channels 
going in (and out of) lead i. It is useful to define 



N„ 



N = Ni+ N2 
mm{Ni,N2) , 



N3 
Na 



\Ni-N2\ 
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The N X N scattering matrix between all three leads can 
be written in block form as 







tl2 


tl3 


s = 


t21 


r22 


t23 




\t31 


t32 


r33 



(3) 



where ru is Ni x Ni and tij is Ni x Nj. A quantity of in- 
terest is the iVi X N2 matrix ti2. Making a singular value 
decomposition (SVD) we have ti2 = UDV where U is 
an Ni dimensional unitary matrix, V is an N2 dimen- 
sional unitary matrix, and D is an A^i x A'2 matrix with 
the (real) singular values c?i > along the Nm diagonal 
elements and zeros elsewhere. More physically relevant 
perhaps is the Hermitian matrix 



T12 =tl2ti2 = V^DtDV 



(4) 



whose N„i nonzero eigenvalues Zi = d^ correspond to the 
squares of the singular values of ti2. It is convenient to 
work in terms of the eigenvalues {zi} rather than the sin- 
gular values {di}. Many physical quantities are given in 
terms of these eigenvalues For example, the current 
out of lead 2 is given by I2 = (712^0 where Go = 2e^//i 
and 



;?12 = Tr[Ti2] - E 



Nrr 



(5) 



Similarly the shot noise of the current at lead 2 can be 
shownjllll to be given by 2eVGns with s = Tr[Ti2(l - 
T12)] = X^ilri ~ ^j)- These physical quantities are 
examples of so-called "linear statistics" whose expec- 
tations can be written as an integral over the eigenvalue 
density p{z) of the nonzero eigenvalues of T12. Other 
linear statistics can also be shown to give physically in- 
teresting quantities related to higher moments of the 
current [3, a £1]. 

We note that a typical element of the matrix S has 
mag nitude (since SS^ = 1), so increasing A^a 



2 



has the effect of decreasing all the elements of ti2. In 
addition to this normalization change, qualitatively we 
should expect that the symmetry (COE, CUE, or CSE) 
constraints on S become less visible in the submatrix ti2 
as -/V3 is increased, since the constraints are more easily 
satisfied by the large number of additional variables in 
the other blocks of S. Thus, we expect that ti2 should 
start to look like an i.i.d. random complex Gaussian ma- 
trix as N3 becomes large. When ti2 is i.i.d. random 
complex Gaussian, the corresponding T12 is known as a 
"complex Wishart" matrix[lll|i|. 

As is frequently the case in random matrix problems, 
it is easiest to make further progress if one focuses on S 
being from the unitary ensemble (broken time reversal 
invariance). In this case we use a method analogous to 
that of Ref. ^3 to determine the joint probability density 
(JPD) of the nonzero eigenvalues {zi} of T12. We start 
by writing the probability density 'P(S) oc (5(8^8 — 1) 
where we have an N'^ dimensional complex S function 
to enforce unitarity. Analogous to 10] we integrate out 
r33,r22,t3i,t32,t2i,t23 to obtain 

P(ti2) K J dti3 j drn S{r\,rn + t\^tr2 + tlgtia - 1) (6) 

with the delta function being complex dimensional, 
and the integrals being Ni{Ni + N^) complex dimen- 
sional. Writing t|2ti2 in terms of its SVD as in Eq. 01 
the Jacobian of this change of variables is given by 
J ~ A({z})^ Yl- z^^ where A({2:}) is the Vandermonde 
determinant 

A(W)- n (^'-^^■) - det[zr']- (7) 

l<i<j<N„^ 



We can then integrate out U and V to obtain 

iTll + t|3ti2 + Z - 



1) (8) 



where z is a diagonal matrix of the Nm nonzero eigen- 
values Zi (with zeros along the remaining diagonals if 
Ni > N2). Performing the remaining integrals yields the 
final result for the JPD 



N„ 



V{{z})<xAi{z}fY[w{z,) 

i=l 

with weight function 



w{z) 



(i-2)^^e(z)e(i-z) 



(9) 



(10) 



where the step function Q(x) is 1 for x > and is zero 
otherwise. In the limit of large N3 (with Na and N,n 
fixed) the factor (1 — z)^^ becomes e^^^^ (^ e~^^ since 
z{Ni + N2) becomes small) as expected for the limiting 
complex Wishart distribution[3,[3|i,0|. 

It is interesting that Eqs. 1^ 1 101 are similar in form to 
the distribution found for the (complex) eigenvalues of 



a (square) truncated unitary matrix found by Ref. |lO| |. 
Here, in contrast, we are looking at the sm^w/ar values of 
the truncated unitary matrix (ti2). Despite the similarity 
of these cases, we will see below that the distribution of 
the singular values of ti2 is very different from that of 
the absolute value of its eigenvalues. 

The Vandermonde determinant can be recast as a 
Slater determinant of a set of wavefunctions {4>i} which 
are orthonormal with respect to the weight function w{z) 



<t>i{z) 



(2z-l) 

,{a,b}. 



(11) 



for i = . . . Nm — 1 where ' (x) is the Jacobi 
polynomial [ni, and Cf = [{2i + N3 + Na + 1)«!(« + 
N3 + iVA)!]/[(i + N^y. [i + TVa)!]- Thus our distribution 
is a realization of the so-called Jacobi Unitary Ensem- 
ble (JUE)[i[0. Interestingly this ensemble occurs in a 
completely different random matrix problem when one 
considers the eigenvalues of A{A + B)~^ where A and 
B are both complex Wishart matrices!^ 0|. In the case 
of the quantum dot with two leads {N^ = 0) it has been 
previously pointed out^j that the eigenvalue distribution 
is a special case of the JUE. Here we have shown that all 
possible JUEs may be realized in three lead dots. 

Using standard techniques 1, 9], once we have our JPD 
in this orthogonal polynomial form, we can write the ker- 
nel K{zi, Z2) as 

N^-l 

E 



K{zi,Z2) = y/ Vj{zi)w{z2) 22 '/'i(^l)<^i(^2) (12) 



i=0 

and this series can be resummed using the Christoffel- 
Darboux formula as usual 1, 9J if desired. All correlation 
functions between the eigenvalues {zi \ can be expressed 
analytically in terms of this kernelll!,l9|, the simplest ex- 
ample of which is the normalized single eigenvalue density 
p{z) — K{z, z)/Nm which is shown in Fig. 1. For com- 
parison, the orthogonal case is also shown (data from 
numerical monte-carlo simulations). In both cases, the 
distribution clearly converges to the complex Wishart 
form0, S P|for large N3. Similar convergence to the 
Wishart form is obtained from numerical simulations for 
the symplectic case (not shown). 

We can also ask about the distribution of 512 — Tr[Ti2] 
which, as mentioned above, corresponds to the conduc- 
tance of our dot. We have 

P(512) = J d{z,} 6ig,2 - Y.^z^) V {{z}) (13) 
= /§/d{^Je-^-(^--i:.-)7'({z}) (14) 

It is now convenient to rewrite the A({zi}) in P({z}) 
in its determinant form (Eq. 01 and expand this deter- 
minant as an antisymmetrized product. We can then 
perform the integrals to obtain 
1 r fih 
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FIG. 1: Top: Scaled eigenvalue density p{zN)/N of T12 = 
tj2ti2- Bottom: Scaled distribution P{gi2N)/N with 312 = 
Tr[Ti2]. In all cases data is shown for A''i = 3,N2 — 2 and 
N3 = 0, 1, 2, 5, 15 (note iV = Ni+N2+Ns). The unitary cases 
are analytic as discussed in the text. The orthogonal cases 
{N3 7^ 0) are numerical and the curves have been smoothed for 
clarity. The orthogonal = case is also analytic as given 
by Refs. Ulli. The dashed curves are the Wishart limit Q 
corresponding to ^ cxd and is the same for both unitary 
and orthogonal cases. The eigenvalue density is only nonzero 
for z < N, and P{g\2N) is only nonzero for 312 < Nm = 2. 



where F is an Nm x Nm matrix with elements 



Fab{k) = dzz 
Jo 



(1 - z)^^ e'^^ (16) 



The integrals of Eq. ^| although tedious can always be 
done, resulting in the form 



where each L„ is an Nm {N — Nm) — 1 degree polynomial. 
The detailed forms of these polynomials can easily be 
evaluated on a case-by-case basis (trivially with Mathe- 
matica). In the limit of large A^3, the (1 — z)^^ factor 
in the integral in Eq. ^|can be replaced by e~^^^ and 
the upper limit of the integral can be extended to 00. We 
then have Fab{k) = (7VA-ha-h6-2)!(iV3-ifc)(i-°-^^^^)- 
The determinant is then detF(fc) ~ (N3 — ik) 



912 



-912^3 



N1N2-I -gr2N 

i/12 ^ I 



suiting in P{gi2) 
which is known as the Rayleigh distribution — the char- 
acteristic distribution for the trace of a complex Wishart 
matrix. 

In Fig. 1, examples of the distribution P[gi2) are 
shown. The results in the unitary case are analytic as de- 
scribed above, whereas the orthogonal data has been gen- 
erated numerically. Convergence to the complex Wishart 
(Rayleigh) limit is again clear in both cases. Similar con- 
vergence to this limiting form is obtained from numerical 
simulations for the symplectic case (not shown). 

As with most random matrix problems, considerable 
simplification occurs in the limit of large matrices. We 
will thus consider the limit of iVi , 7V2 , -^3 — * 00 while 
keeping the ratios Ni/N fixed. We proceed by consider- 
ing Eq. 1^1 as a classical one-dimensional gas with Boltz- 
mann weight given by "Pdz}). Thus, we have a one par- 
ticle potential given by e"^*-^-* — w{z), and two parti- 
cle interaction e~^^^''^^^ = {zi — ZjY . Using standard 
methods 0,|i,E3, we derive the saddle point equations 



We then introduce the Green's function 



(20) 



(21) 



with p{z) — TT^^ IniG'(z + iO^) being the normalized sin- 
gle eigenvalue density. It is easy to verify P, 0,0] that to 
order 1/A^ we have G{z)'^ = ^ V'{zi)/{z — Zi) which we 
can solve analytically resulting in the eigenvalue density 



Fab{k) = Qab{l/k) + e'''Rab{l/k) (17) 

where Qij and Rij are simple polynomials of degree a + 
6—1-1- A^A + -^3, and Fij (0) remains finite. We then have 



p(z) 



J^^{z-z^){z+~z) 

Nm 



for z_ < z < z+ and p{z) 



2ttz{1 - z) 
otherwise, where 



(22) 



det[F(fc)] = ^i?„(l/fc)e^^-" 



(18) 



n=0 



where the D^s are polynomials of order Nm{N — Nm)- 
The Fourier transform in Ea. llSl can easily be carried out 
analytically resulting in P{g) being of the form 



P{9i2) = E ■^"(512)0(512 - n) 



(19) 



z±=N-^ (^^Ni{N^N2) ± ^N2{N-Ni) 

(23) 

We have also derived this simple form by using the free- 
probability Ly| method, which we will present in a sepa- 
rate paper |l4). In fact, we also show in that work that 
this form is correct in the large N limit independent of 
whether the original scattering matrix S is from the COE, 
CUE, or CSE. This result should not be surprising, as it 
is well known^ 0] that for many quantities of physical 
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FIG. 2: Eigenvalue density p{z) of T12 = t|2ti2 for iVi/iV = 
1/6, Ar2/A^ = 3/6 and N3/N = 2/6 for iVi = 1,2,3,5 
(A'^ = 6,12,18,30). The unitary case is calculated analyti- 
cally as described in the text. The orthogonal case is calcu- 
lated numerically. For A'^i = 1,2 the orthogonal case numer- 
ical data has been smoothed for clarity. For A'^i = 3 (circles) 
and Ni = 5 (pluses), smoothing is found to destroy the oscil- 
lations, so the raw histogrammed data is shown. The dashed 
line is the large A'^ limit (Eq. 1221 to which both the unitary 
and orthogonal cases converge. 

interest (such as the average conductance, etc), the en- 
sembles differ from each other only at order unity whereas 
leading terms will be order N. These small differences 
are known as the weak localization corrections. 

In Fig. 2 we show how the eigenvalue distributions con- 
verge to the large N limit, holding Ni/N fixed. Analytic 
results are shown for the unitary case, and numerical re- 
sults are shown for the orthogonal case. In both cases 
we see convergence to the limiting form of Eq. |23 as ex- 
pected. Similar convergence to this limit is seen numeri- 
cally for the symplectic case (not shown). 

As mentioned above, the problem of calculating the 
eigenvalue distribution of T12, or equivalently the singu- 
lar value distribution of ti2 appears similar to the prob- 
lem of calculating the complex eigenvalue distribution of 
a (square) truncated random unitary matrix derived in 
To see the difference between these two problems, 
let us examine the case of Ni/N — N2/N = n < 1/2 
in the large N limit for which both calculations are 
applicable. From our results, we find the squared sin- 
gular values density decreases monotonically to zero at 
z — z-^- — Anil — n). In the work of [inj |. however, it is 
found that the density of the norm squared of the eigen- 
value distribution is a function that increases monotoni- 
cally, then drops to zero discontinuously at z — n. 

One might attempt to guess at the general form of the 
JPD, analogous to Eqs. ElandUm for the COE and CSE. 
Indeed, in the case oi — e. simple generalization 
exists 0, ^3 which involves changing the exponent of z 



in the weight ^(z) and changing the exponent of the 
Vandermonde determinant in Eq. |3 to /3 = 1 for the 
COE and to /3 = 4 for the CSE case. However, for the 
-^^3 7^ problem no such simple form can apply, since 
for large N3 we always expect that T12 will be Wishart 
(ie, ti2 is simply an i.i.d. complex Gaussian matrix) so 
that there will always be quadratic (/? — 2) eigenvalue 
repulsion independent of whether S is from the COE, 
CUE, or CSE. We thus leave the full solution of the JPD 
for the COE and CSE as an open problem. 

The authors are greatly indebted to A. Lamacraft and 
B. Simons for pointing us to this crossover as an inter- 
esting problem to examine. 
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